Extrinsic levels, diffusion, and unusual incorporation mechanism of lithium in GaN 
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Results of a first-principles study of the Li impurity in GaN are presented. We find Li is 
a channel interstitial, with an onset for diffusion at T~ 600 K. Above this temperature, Li can 
transform to a Ga-substitutional acceptor by exothermic recombination with Ga vacancies. This 
process implies capture of at least one electron; therefore Li acts as an electron sink. Lioa is stable 
again interstitialcy, and has a shallow first ionization levels of 0.16 eV, and second ionization at 
0.63 eV. Lattice locations and their temperature dependence are in close agreement with recent 
experiments. 

PACS numbers : 71.25.Eq, 61.72.Vv, 61.72.Ss 
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Doping IH-V nitrides p— type is a conceptually intrigu- 
ing and practically relevant problem. The current practi- 
cal recipe for GaN p— doping is based on the Mg acceptor, 
but the search for alternative dopants is still ongoing. 
Useful informations on the mechanism involved in Mg 
incorporation Jd and on possible alternative acceptors ,□ 
has come from theoretical investigations^ While Group- 
II dopants have been treated in detail,Era little atten- 
tion has been devoted so far to potential Group-I double 
acceptors Jj 

In this Communication we discuss the physics of Li as 
an impurity in GaN, building on direct calculations of 
energetics, equilibrium geometries, and diffusion paths. 
Li is found to be preferentially incorporated in GaN as 
an interstitial. Since its chemical incorporation is costly, 
incorporation should be effected by ionic implantation. 
Diffusion of Li in GaN is activated, with a threshold 
temperature of about 600 K, and occurs in three dimen- 
sions with a slight preference for motion within the basal 
plane. Li recombination with Ga vacancies is a highly 
exothermic, barrier-less reaction. Hence, if diffusion is 
activated, Li is efficiently incorporated substitutionally 
into any available Ga vacancy. In the recombination pro- 
cess, capture of one or two electrons occurs for almost all 
values of the Fermi level Ep. The substitutional double 
acceptor Liga thus realized is highly stable again inter- 
stitialcy (as the Ga vacancy-Li interstitial pair is highly 
disfavored with respect to Lica) and, most importantly, 
it exhibits a thermal first ionization energy of 0.16 eV, 
smaller than the ionization energy of Mgca, and a second 
ionization level at 0.63 eV. Because of the strong depen- 
dence on precursor defects (the Ga vacancy), the role of 
Li as an acceptor or compensator will be significant only 
if its precursors are the dominant native defects of the 
system. Of course, this limits its practical relevance for 
GaN doping. 

We used local density-functional theoryO ultrasoft- 
pseudopotentialQ plane-wave calculations of total ener- 
gies and forces in GaN wurtzite supercells (typically en- 
compassing 32 atoms), and with plane wave cutoff of 
25 Ry, to predict from first principles the total energies 



and equilibrium geometries relevant to our problem. All 
atomic coordinates are fully relaxed until residual forces 
are below 0.01 eV/A. Further technicalities and the for- 
malism to calculate defect formation energies have been 
reported previously in detail (see e.g. Refs. ^and ||). One 
point worth mentioning here is the solubility limit being 
metallic bulk Li, rather than LisN as could be expected. 

The formation energies of relevant defects and impu- 
rity configurations are summarized for reference in Fig- 
ure |l|. We begin by examining the substitutional Lica- 
In this configuration Li behaves as a double acceptor: we 
determined directly its first and second thermal ioniza- 
tion energies, that are 0|J6 eV and 0.63 eV, respectively, 
above the valence bandE This gives us additional moti- 
vation for this study, since the first ionization energy is 
comparable than that of Mgca! the second ionizatioruia. 
fourfold the first as expected from hydrogenic models .H'cl 

Since the chemical formation energy of Lica is rather 
high (3.83 eV in the neutral state in N-rich conditions), 
plain chemical solubility does not seem an easily viable 
incorporation pathway for substitutional Li in rGaN in 
moderately p— type conditions. Ion implantationtd could 
possibly surrogate chemical reaction. It is then natural to 
investigate precursors to substitutional incorporation, se- 
lected among native defects and Li complexes therewith. 
Specifically, we studied the Ga vacancy and a selection 
of Li interstitials and Vca-Liint pairs. We also mapped 
out the diffusion path and estimated the activation bar- 
riers for Liint motion along high symmetry directions in 
undefected wurtzite GaN, and towards a neighboring Ga 
vacancy in cation-vacancy defected GaN. 

The results for VcarP-nd Vn are quite compatible with 
previous calculations.u The cation vacancy, in particu- 
lar, is a deep triple acceptor, with formation energy go- 
ing from fairly high in p— type conditions to very low in 
n— type. Local relaxations are in agreement with previ- 
ously reported data. 

As for Liint in undefected GaN, we looked for stable 
sites in the open channels of the wurtzite structure and 
in the adjacent trigonal cages. In all configurations the 
stable charge state is +1. We found that the marked 
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preference [for antibonding or bond-center configurations 
exhibitedBlij by hydrogen in GaN is not shared at all by 
Li. None of the various sites suggested by a naive analogy 
with H is found to be stable. The lowest-energy configu- 
ration of Liint is at the center of the open wurtzite chan- 
nel, and has C^y symmetry. The site can be viewed as 
delimited at the top and bottom by three-fold asymmet- 
ric boat-shaped Ga-N rings, and laterally by symmetric 
boat-shaped three- fold rings. The coarse-grained picture 
is that Li simply tries to maximize its distance from all 
other atoms, but on closer inspection a slight preference 
is revealed towards N atoms. In particular, Liint sits at 
3.61 bohr from the N atoms protruding from the upper 
ring, and at 4.33 bohr from the Ga atoms in the lower 
ring. The equilibrium sites along the c axis ara-0.23 c and 
0.73 c, which compare well with the measuredll3 values of 
0.25 c and 0.75 c [the implicit convention here is that Ga 
atoms have z-coordinates zero and 0.5 c, and N atoms 
have 0.376 c and 0.876 c ]. Inside the trigonal cage, we 
only find a marginally stable site near the cage center, 
at about 1.3 eV higher than the channel site, and with a 
confinement barrier of only 0.1 eV. Again, Li prefers to 
sit closer to N atoms (Li-N distance ~ 3.18 bohr) than 
to Ga atoms (Li-Ga distance ~ 3.59 bohr). 

These results can be largely understood with a sim- 
ple Madelung-like electrostatic model, whereby we cal- 
culate the electrostatic energy of channel-interstitial Li 
as a function of its coordinate along the c axis, assuming 
it interacts with point-charge ions placed on a wurtzite 
latticej— I N and Ga are assigned their dynamical Born 
chargdlj of ± 2.7, while Li is assumed to have a dy- 
namical charge of 0.9. The equilibrium site is found at 
0.31 c, close to 0.23 c and 0.25 c given by direct calcu- 
lation and experiment, respectively; the barrier for dis- 
placement along c is 3 eV, which compares well with 2.7 
eV from direct calculation for the unrelaxed lattice. This 
shows that basically Li sits at the site dictated by the 
minimization of the electrostatic interactions with the 
ionic lattice. The sizable deviation in the equilibrium co- 
ordinate is due to the large size of the N ions, assumed 
to be pointlike in the model. 

Starting from these sites, we investigated the diffu- 
sion paths for Li in GaN. The barrier for diffusion be- 
tween equilibrium sites along the (0001) direction is 1.55 
eV, and the vibrational frequency at the equilibrium site 
along (0001) is 13 THz. When the diffusion path from 
channel to channel across a cage is considered, the max- 
imum barrier is found to occur near the cage center, at 
1.4 eV above the channel site. At temperatures at which 
this diffusion path is activated, the previously mentioned 
marginally stable pocket near the cage center is of course 
irrelevant, given its confinement barrier of only 0.1 eV. 
The vibrational frequency at the channel equilibrium site 
along the cross-cage diffusion direction is 9 THz. The ef- 
fective activation temperature for Li diffusion can be es- 
timated from the above data to be about 600 K. At the 
activation temperature, cross-cage motion, i.e. motion 
approximately within the basal plane, is more probable 



than motion along the c axis by a factor of 10 on account 
of its lower barrier. 

We now consider the propensity of Li to recombine 
into Ga vacancies, which we assume henceforth to be al- 
ready present in the ^stem. In Figure |2| we compare a 
distant Li^-Voa pairEj with substitutional Lioa; the lat- 
ter is always vastly favored energetically (by over 5 eV 
for a semi- insulating crystal, i.e. midgap Fermi level), 
so Li- vacancy recombination is highly favorable. This 
recombination process involves several steps. First, Li 
diffuses across the crystal towards the neighborhood of 
a vacant Ga site; as discussed previously, this requires 
a temperature above about 600 K. Second, Li recom- 
bines with the vacancy: upon direct calculation, we find 
that Li drops into Vca from the wurtzite channel with- 
out any activation barrier, irrespective of the vacancy 
charge state. Third, the substitutional Li thus produced 
assumes a charge state consistent with the Fermi energy. 

As to the latter step we note that, on account of the 
charge state of the precursor Ga vacancy and depend- 
ing on the Fermi level, the metastable Vca-Liint distant 
pair exists in charge states from -1-1 to -2. Upon recom- 
bination with the vacancy, Lica may therefore need to 
capture electrons to reach its equilibrium charge state (0 
to -2, but in fact equal to the latter for most of the ac- 
cessible Ep range). The level diagram in Fig. || shows 
indeed that either one or two electrons will be captured 
over practically all the Fermi level range. For example, in 
process A the pair is singly negative (Li recombines with 
Vq^ ) 7 hence Li captures one electron to reach its equilib- 
rium charge state of -2. In process B, Li recombines with 
^Ga (neutral pair) and capture of two electrons ensues. 
Processes in p-type conditions, such as C in Fig. |^, are 
complicated by the acceptor levels of Lica, but all imply 
the capture of either one or two electrons. 

In short, the reason is that all charge states of the Li- 
vacancy pair overlap in energy only with more negative 
charge states of Lica- The only exception is the doubly 
negative Li-vacancy pair (Li recombination with triply 
negative Ga vacancies): the latter, however, only occurs 
in extreme n-type conditions, where in fact the direct 
chemical incorporation of Li becomes favorable. For all 
other Fermi levels, Li acts as an electron sink, i.e. as an 
effective acceptor, or more properly, compensator. Fig. 
H also shows that electron capture is more likely at lower 
Fermi levels, i.e. towards p-type conditions: in partic- 
ular, if Lioa were always -2, the number of captured 
electrons would change from through 3 as Ep moves 
downwards across the gap. This is quite clearly the op- 
posite of what happens to normal "chemical" acceptors. 
Unfortunately, the mechanism is limited by the availabil- 
ity of Ga vacancies as precursors; any Li not paired to a 
Ga vacancy will remain interstitial. 

Let us summarize the picture we arrived at. Isolated 
Li in GaN sits in the wurtzite channels. It will be very 
strongly attracted towards any neighboring Ga vacancy, 
and it will recombine with it to give rise to substitutional 
Lica- To reach a vacancy, Li must diffuse into GaN; dif- 
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fusion sets in at approximately 600 K (with a preference 
for the basal plane near the onset). Once in the neigh- 
borhood of the vacancy, Li recombines with the vacancy 
without further energetic barrier hindering the process. 
Given the gain upon recombination of Li into the va- 
cancy, the reverse process Lica Lijnt + Vca is com- 
petely ruled out, i.e. Lica is stable again interstitialcy. 
At the end of the recombination transient, the Lioa con- 
centration will equal that of the pre-existing Ga vacancies 
(which will have vanished), and the electron density will 
have been reduced by one to two times the smaller of 
the two concentrations of Lijnt and Voa, depending on 
the initial Fermi level. Clearly, compensation of n-type 
carriers will effectively occur if Li is incorporated in GaN 
in its ionized state, as in ion implantation. Any Li 
population in excess of the Ga vacancy concentration will 
remain interstitial. 

We close our investigation making contact with the 
only piece of experimental evidence on Li in GaN, 
a study of [decay emission from implanted radioactive 
Li isotopes. liEl At room temperature, channeling 7— ray 
emission is observed, showing that Li is indeed intersti- 
tial; as mentioned earlier on in this paper, our calculated 
locations agree closely with those deduced from angle- 
dependent channeling. After annealing at 700 K, the 7 
emission from Li decay is shadowed, indicating that Li 
has either moved into a trigonal cage, or to a substitu- 
tional site. On the basis of our results, the natural expla- 
nation is Li-Ga vacancy recombination. While shadowing 
may occur also if Li would remain confined in a trigonal 
cage, our calculated diffusion barriers show that this is 
not possible, as the cage center is not a stable site at the 
relevant temperature. Our estimate of 600 K for the on- 
set of diffusion (and consequent recombination with Ga 
vacancies) is in rather good agreement with experiment, 
given our neglect of entropic contributions to the free 
energy diffusion barrier. 

In summary, we have shown by first-principles calcula- 
tions that Li can be easily incorporated as an interstitial 
in GaN, and that it will efficiently recombine with cation 
vacancies to yield Liq^. The resulting substitutional dou- 
ble acceptor is predicted to have first and second thermal 
ionization energies at 0.16 eV and 0.63 eV above the va- 
lence band, energies which are comparable or smaller to 
those of acceptors commonly used in GaN. Due to lim- 
itations due to the availability of its precursor Vca, Li 
is expected to function only as compensator of n-type 
carriers. 
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